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Abstract 

We calculate the flux of ultra-high energy protons due to the process of “cusp evaporation 
from cosmic string loops. For the ‘standard 'value of the dimensionless cosmic string param- 
eter t = Qfx ~ 10 -6 , the flux is several orders of magnitude below the observed cosmic-ray 
flux of ultra-high energy protons. However, the flux at any energy initially increases as 
the value of e is decreased. This at first suggests that there may be a lower limit on the 
value of £ , which would imply a lower limit on the temperature of a cosnuc-string-forming 
phase transition in the early universe. However, our calculation shows that this is not the 
case — the particle flux at any energy reaches its highest value at e ~ 10 15 and it then 
decreases for further decrease of the value of e. This is due to the fact that for too small 
values of e(< 10~ 15 ) t the energy-loss of the loops through the cusp evaporation process 
itself(rather than gravitational energy-loss of the loops) becomes the dominant factor that 
controls the behavior of the number-density of the loops at the relevant times of emission 
of the particles. The highest flux at any energy remains at least four orders of magnitude 
below the observed flux. There is thus no lower limit on e. 

( V A " A — C 1 — 1 c ’ 9 ? h ) • T r t ‘ ' U n J- 1 1’*+ "V 

UL A-MliJH P|\r ; <GY COS'UC f«Y, (O’ ic^ so 

uni v. ) I f~: i* t . L ^ J 

Unci is 

Operated by Unlvereitie* Reeearch Association Inc. under contract with the United States Department of Energy 



value of € - — an assumption, which, as we have mentioned above and shall discuss below, 
is not valid. We will report the explicit calculations for the case of neutrinos elsewhere, 
but from the results of Ref. 4 and the discussions given below, it already appears that the 
use of the correct formulas for the loop number densities would also eliminate the lower 
bounds on e found in Ref. 4. 

In Sec. II, we briefly describe the process of cusp evaporation from CS loops and 
estimate the number of primary particles emitted from the string per unit time. The 
UHE proton injection spectrum, resulting from the decay of the primary particles and the 
subsequent hadronization of the decay products, is estimated in Sec. Ill by using a suitable 
hadronic jet fragmentation distribution function. A general expression for the predicted 
flux in the present epoch is written down in Sec. IV. In Sec. V, we briefly discuss the 
main processes by which UHE protons lose energy during their propagation through the 
cosmic medium, and discuss how the effective maximum possible redshift of injection is 
determined for a given value of the energy of the proton in the present epoch. The CS loop 
length distribution function required for our calculation is obtained in Sec. VI. The main 
calculation of the flux is described in Sec. VII, and the results, discussions and conclusions 
are presented in Sec. VIII. 

Except where otherwise stated, we use natural units, ft = c = 1, so that \fG = = 

tpi, where Mpi is the Planck mass and tpi is the Planck time. The Hubble constant is 
Hq = 100. ft Km.s~ l .Mpc -1 , and we use ft = 0.75. Also t eq is the time of equal matter- 
and radiation energy density, z eq is the corresponding redshift and <o is the present age of 
the universe. We assume a fio = 1 universe. 

II. CUSP EVAPORATION 

A non-self-intersecting, 5 ’ 6 freely oscillating CS loop has one or more points which momen- 
tarily achieve the speed of light once during every oscillation period. These points called 
“cusps” appear 6 if the motion of the loop is described by the Nambu action, which is valid 
for infinitely thin strings. In reality, CS have a finite width, and so the Nambu action 
is, strictly speaking, not valid for CS and true cusps may not form. Nevertheless, “near 
cusp” points are likely to occur where the string moves with very high Lorentz factor. 
At a cusp, two string segments overlap, and it has been pointed out 3 that interactions 
of the underlying fields lead to ‘evaporation ’of the overlapped region whereby the energy 
contained in the overlapped region of the loop is released in the form of particles, thus 
smoothing out the cusp. New cusps continue to form and evaporate during each period 
of oscillation of the loop. The length of the cusp region of the loop can be estimated 3 as 
^cusp ~ where L is the total length of the loop and u> ~ /z -1 / 2 is the width of 

the string. (The length L of the string is defined such that pL is equal to the total energy 
of the string). The energy released due to cusp evaporation will be in the form of bursts 
with time-scale A tburst ~ ^ cusp • The period of oscillation, T osc , for a loop of length L 

l/3 

is 5 y. Thus, Atburst /T osc -(«'• «i- Thus the rate of energy released due to cusp 
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reproduces the particle multiplicity growth as seen in GeV-TeV jets in colliders. This 
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gives, 


f -IS 


(4) 


where x = E/Ej e t — e~jh < 1, E being the energy of a hadron in the jet. A small 
fraction(~ 3%) (Ref. 11) of the hadrons in the jet will be nucleons and antinucleons which 
ultimately end up as protons and antiprotons. Observationally, since the primary particles 
at the high energies involved here are not detected directly, one cannot distinguish between 
protons and antiprotons. We shall, therefore, in the following, collectively refer to them 
simply as protons. Let denote the injection spectrum of the protons, i.e., the 

number density of injected protons per unit energy interval at an injection energy E, per 
unit time at an injection time t, due to cusp evaporation from all CS loops. Then using 
eqs.(4) and (2) we get, 


*(£,,*,) a 2 x 0.03 [ dL^(L,ti)L~3, (5) 

where x — 3 E{jE\ = 3 E{f~ x e~^Mp] and j y(L,ti ) is the CS loop length distribution 

dJu 

function, i.e., dn(L,ti) is the number density of CS loops with lengths in the interval 
[L, L + dL] at the time t,. The factor of 2 in eq.(5) takes care of the fact that we have 
assumed two quarks in the decay products of each X and each quark produces one hadronic 
jet. Thus eq.(4) yields an injection spectrum oc E~ 3 ^ 2 for Ei « Ex- 


IV. GENERAL EXPRESSION FOR THE FLUX 

Let j(Eo ) denote the number of protons per unit energy interval at energy Eq in the 
present epoch(to) crossing per unit area per unit solid angle per unit time due to the 
source $(Ej,<i). Then, assuming an isotropic distribution of the CS loops in an Einstein- 
deSitter “flat” (Do = 1) imi verse, we get 

OO 

XEo) = ±J*™ 3 (t i ydr[(l+z i )-'i(E i ,t i )] (g)^ W( l (o)r2 . (6) 

o 0 

where t, is the injection time, z, is the corresponding redshift, Ei = E t (Eo.t,) is the energy 
at the time of injection tj, a(t) is the scale-factor of the universe, and r is the comoving 
radial coordinate of the source. The factor (1 + z,) _1 = a(ti)/a(to) in eq.(6) is due to 
the cosmological “redshift” of the frequency of emission. 14 Now for a fio = 1 universe, 
r = c f*° dt/a(t) (assuming that the particles are ultrarelativistic, so that they travel 
almost with the speed of light, c), so that a(ti)dr = — cdt{ . Furthermore, t , > t eq (in fact, 
as we shall see below, for all values of energy Eq, all injection times f, satisfy ti >> t eq ) 
so that (1 + Zi)~ l = a(ti)/a(to) = (ti/to) 2 ^ 3 , giving a(t{)dr = — cdti = §cto(l + z t )~ 5 / 2 dzi. 
Putting all these together, eq.(6) becomes 

j(E o) = jT <{*(! + ZiT"' 2 ( ^|r ~ ) e (7) 
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cosmic-rays, while not entirely devoid of controversies, do seem to indicate 19 the existence 
of a cutoff as predicted. 

Now, given the full knowledge of the energy-loss functions /?o , P air{E) and Po,pion{E), 
one can solve eq. (11) numerically to find the energy E, of a proton at any injection 
redshift z t corresponding to a given value of its energy in the present epoch(£’ 0 ). One can 
then evaluate the injection spectrum ) using eq. (5)(with a given CS loop length 

distribution function, see Sec. VI) and obtain the flux by evaluating the z, -integral in 
eq. (7). The full numerical calculation according to this procedure is described in Ref. 20 
in the context of another particle production process involving CS. Here we undertake 
an approximate calculation which essentially yields the same result, but it allows us to 
avoid the full numerical solution of eq. (11). The approximation is based on the use of the 
arg um ents that lead to the prediction of the Greisen-Zatsepin-Kuz’min 18 cutoff mentioned 
above. 

To see this, let us consider the energy-range 5 x 10 18 eV ^ f?o ^ 6 x 10 19 eV r , in which, 
as mentioned above, /3 0 , pair is dominant over do.p.on and the former is weakly energy- 
dependent remaining roughly constant at /3q ss 2.13 x 10 -1 °yr _1 . In this case, as long as 
(1 + Zi)E{ < 6 x 10 19 eV, eq. (11) has the analytic solution, namely, 


Ei(zi) = £o(l + Zi) exp 


3 Ho 



3 

+ Zi )2 



for 5 x 10 18 eV ^£ 0 ^6x 10 19 eV. 


( 12 ) 

Thus in the above energy-range, if we consider a proton at energy Eo today, its energy E, 
at any injection redshift z, rises exponentially with z t . If for any given value of Eo , we 
define the injection redshift z , ima x such that 


(1 + )Ei(zi 


— Z i, max') 


E 0 ) ~ 6 x 10 19 eV, 


(13) 


then for z* > z, mar , the proton would be in the photopion energy-loss regime. In this 
regime the energy-loss itself rises sharply (roughly 17 exponentially) with energy and so the 
energy E{ of the proton at the injection redshifts Zi > Zi tTnax (Eo ) rises evan faster 21 with 
increasing values of z t . As a result, the rapid fall of the injection spectrum ( $(E t , z t ) (which 
goes as ~ E~ 3 ^ 2 ) with increasing value of z t dominates over the power- law rise of $ with z, 
coming from the fact that the number- density of the CS loops increases with redshift(see 
Sec. VI-VII). This in fact ensures that the . 2 ,-integral in eq. (7) converges fast. In other 
words, for a given value of E 0 , the contribution to the flux j(E 0 ) of eq. (7) from injection 
redshifts Zi > ~z limax (Eo) are negligible compared to those from the injection redshifts 
Zi < Zi tTnax (Eo). The quantity z, mai defined by eq. (13) can, therefore, be taken as an 
effective cutoff for the integral in eq. (7). Actually, since the maximum energy of a particle 
in our case cannot exceed 1 Ex , the cutoff redshift should be determined from the condition 


(1 + Zi jmax )Ej(z, — Zi <max , Eo) — min ^£lx(l + z ljTnax ), 6 x 10 eV^ , 


(14) 
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(ii) for t > t*, 


if Tet < L < at, 


^ L dL ' 

dn(L,t) - j ^(p^j-4 ('ffffiMy dL, if LZ”(t) <L< Tef, <19> 

,0, ° if L < L^in- 

In deriving eq.(18)-(19) we have assumed that the loops survive with their lengths essen- 
tially unchanged till the end of their lifetime at which they instantaneously disappear. 

VII. CALCULATION OF THE FLUX 

We are now ready to evaluate the flux from eq. (7). First let us define Ii = 

C /3 f dL—(L,ti)L ~ 1 ^ 3 . Using (18) and (19) together with the appropriate forms for the 
scale factor of the universe in the matter- and radiation dominated epochs, these L-integrals 
are easily carried out. After some algebra, and expressing ti in terms of Z{ by the relation 
ti = to(l + Zi)~ 3/2 , we get 
(i) for ti < t*, 

f K l(l + z i ) 9 ^ 2 — k 2(1 + z i) 5 i f° r (1 + z i) < Z\ , 


Il = 


\ k 3 (1 + 2 ,) 81/16 + k 4 (1 + z , ) 3 - /c 5 (l -I- 2 i) 5 , for (1 + Zi) > Zi, 


( 20 ) 


where Z x = 7 c^ 3 e x ^a 8 / 9 ( tpi/t 0 ) 2 ^ (1 + z eq ) 4 / 3 , and 


«i = X e 1/6 

« 2 = \pa~ x/z , 

4 


tpi 

to 


- 1/3 


/c 4 = -pa - 1 ' 3 ( M 
44 K \t 0 J 


-t/3 


( 21 ) 


« 3 

and 


= A^ a 3/2 7 -ll/8 e ll/48 / ^iV /2 ( f _Pi\ 

11 \ to J \ to J 


- 11/24 *5 — « 2 ; 


(ii) for ti > t+, 

■ __ f -Ai(l + Zi) 5 — ^ 2(1 + Zi) n / A , for(l + Zi) < Z 2 , 

1 £l(l + z i ) 23 ^ 4 — £ 2(1 + z i ) 6 + £ 3(1 + ^i) 3 — £ 4(1 + z i) 5 > f° r (1 + z i ) > Z 2 , 

( 22 ) 


where Z 2 — (re/o) 2 / 3 (l -f z eg ), and 
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One might think, by looking at eq. (24), that decreasing the value of f{i\e.,Ex; see 
eq. (3)) may give a higher value of the flux at any given energy. This is true as long as 
the values of / and e axe such that E x > 6 x 10 19 eF. However, if / is made too small, 
eventually one gets Ex < 6 x 10 19 eV, in which case eq. (14) yields a smaller value of 
z<,mai(than what one would get for the case E x > 6 x 10 19 eV; see Fig. l)resulting in a 
smaller flux. Moreover, for too small values of / and e one gets Ex < Eq, in which case, 
obviously, no particles of the given energy can be produced in the first place. Explicit 
calculation shows that the peak flux always remains below the value obtained with / = 1 
and e ~ 10 -15 . 

Note also that in all the above calculations we have assumed that the cusp evaporation 
process occurs at the maximum efficiency(7 c = 1). If 7 c << 1, then all the above values 
of the fluxes will be correspondingly lower. 

Now, consider the case of neutrinos. First note that for e < 5.43 x 10 -16 , eq. (17) 
gives <* > to ~ 2.67 x 10 17 sec (for Qo = 1, h = 0.75). In this case, obviously, all the 
injection times f, satisfy f, < t* irrespective of whether one is considering neutrinos or 
protons. Eqs. (20), (21) then imply that the values of the flux at all energies will decrease 
with further decrease of the value of e for e < 5.43 x 10 -16 . So the lower limit, e > 10 -17 , 
found in Ref. 4 will probably disappear when the correct form of the loop LDF is used. 
Similarly, for the case e = 10 -15 , we have t* ~ 1.7 x 10 16 sec, and with E x = 10 15 GeV and 
for Eo = 10 19 eV, say, we have for neutrinos, 4,25 1 + Zi,max = Ex/Eo — 10 5 , implying that 
the earliest possible time of injection(#i >m i„) satisfies ti iTntn « t *. So, the contribution to 
the present-day flux from the tj’s in the range t t rnin < f, < <*, when calculated by using 
the loop LDF as determined by cusp evaporation itself(eqs. 20-21), will give a lower value 
of the flux(at the given energy) than what is obtained in Ref. 4. Further reduction of the 
value of e will then reduce the flux further. Thus the lower bound, e > 10 15 , found in 
Ref. 4 will also, it seems, disappear, unless the values of some other parameters(e. g., (3) 
are significantly different from their currently favored values. 

In summary, we have estimated the UHE proton flux resulting from CS cusp evap- 
oration process and found that the flux at all energies remains below the observed flux, 
and that there is no lower limit on the temperature of a CS-forming phase transition in 
the early universe as far as high-energy particle production from CS cusp evaporation is 
concerned. 
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